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Abstract. We investigate static and dynamical ground-state properties of the two- 
impurity Anderson model at half filling in the limit of vanishing impurity separation 
using the dynamical density-matrix renormalization group method. In the weak- 
coupling regime, we find a quantum phase transition as function of inter-impurity 
hopping driven by the charge degrees of freedom. For large values of the local 
Coulomb repulsion, the transition is driven instead by a competition between local 
and non-local magnetic correlations. We find evidence that, in contrast to the usual 
phenomenological picture, it seems to be the bare effective exchange interactions which 
trigger the observed transition. 
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1. Introduction 

Although forty years have passed since the discovery of the Kondo effect, it is still 
one of the most interesting topics in condensed matter physics; it lies at the heart 
of understanding strongly correlated electron systems The Kondo effect, which 
leads to the quenching of an impurity spin, forms the basis of the physics of a single 
magnetic impurity embedded in a metal. However, systems with more than one impurity 
are considerably more complicated and present additional difficulties in a theoretical 
investigation. In particular, there are two effects which compete against each other 
in multiple-impurity systems: the Kondo effect and the Runderman-Kittel-Kasuya- 
Yosida (RKKY) interaction. The RKKY exchange favors the formation of non-local 
magnetic correlations; the Kondo effect, on the other hand, is based on purely local 
magnetic correlations. The competition between Kondo effect and RKKY interaction 
is thought to be the key mechanism to understanding the magnetic properties of the 
heavy fermion materials |2]. The simplest systems in which to study this competition 
are two-impurity models. Recently, such models have also attracted much attention in 
the context of double quantum dots [2], which can be viewed as a direct experimental 
realization of the two-impurity Kondo model. 

Theoretically, the two-impurity Anderson model (TIAM) [113011111113110111111121 
Cni El eg HE] and the two-impurity Kondo model (TIKM) [ni[THllinil2nilZIll221l2Sll211 
l2Sll2nil2Ill2Hll2niEniEIlE21E3IMI have been extensively studied with various methods; 
nevertheless, their physical properties at low temperature are not yet well understood. 
In particular, the situation is far from clear concerning dynamical properties because 
only a few methods are able to reliably calculate the dynamical properties of such 
models due to their complexity. So far, the spectral density has been calculated using 
perturbation theory (PT) P?lH l4j and the numerical renormalization group (NRG) |22j . 
but the results are not fully satisfactory. The PT provides an explicatory and accurate 
picture of quantum impurity dynamics only in certain limiting cases. While the NRG 
can determine the low-energy dynamics of quantum impurity models almost exactly, it 
is less precise at high energy. 

Recently, the dynamical density-matrix renormalization group (DDMRG) 
method |3S1 was applied to the single-impurity Anderson model (SIAM) and it was 
shown that the method can calculate the impurity spectral density with good resolution 
for all frequencies and coupling strengths [HEl EZj- This method can be extended to 
investigate the dynamics of a two-impurity problem without difficulty. Here we study 
the spectral density of the TIAM using the DDMRG method. Since the parameter space 
of the TIAM is rather large, for simplicity we focus here on the limit of small inter- 
impurity distance. This simplification does not change the substance of the problem 
and is quite likely relevant for typical experimental situations, e.g., clusters of magnetic 
atoms on metal surfaces or multiple quantum-dot systems. 

The aim of this paper is to demonstrate the efficiency of the DDMRG for the two- 
impurity system and to discuss the dynamical properties of the TIAM in the limit of zero 
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impurity distance as the first step of a more general DDMRG study. Tlie organization of 
this paper is as follows. In section 2, the model Hamiltonian for the TIAM is introduced. 
In section 3, we transform the Hamiltonian for an efficient treatment with the (D)DMRG 
and define even- and odd-parity orbitals of coupled impurities. In section 4, we show 
the static and the dynamical properties calculated with the (D)DMRG method. The 
conclusion and discussion follow in section 5. 

2. Model 

The Hamiltonian for two impurities placed at Rj {i = 1, 2) is written as 

licr ikcr 

+ UJ2 « - - + ^12 J2idlJ2. + h.c) 

i=l,2 (T 

where d\„ {d^^) creates (annihilates) an electron with spin a =t,i in a local level (the 
impurity site i), hf^ = dl^d^^and fl (/^) creates (annihilates) an electron with spin a 
in an eigenstate of the (noninteracting) host band with dispersion ^k- The sum over 
k runs over all states of the host band. The hybridization between the local impurity 
state and the delocalized band state k is given by the positive couplings Vk. Electrons 
in the local level are subject to a Coulomb repulsion U . In this paper, the energy level 
of impurity sites is set by = — f//2. Under this assumption, we can map our model to 
the two- impurity Kondo model in the strong-coupling limit. We also set Ri — R2 = 
and replace k ^ for simplicity. Thus, our model depends on the parameter U and 
the hybridization function 

A(^) = 7r^ 1141^5(0; -£fc)>0. (1) 

k 

For a symmetric hybridization function, A(u;) = A(— u;), the TIAM is particle-hole 
symmetric for ti2 = 0. 

Since, for the time being, we are interested in understanding the qualitative aspects 
of the model, it is convenient to choose a fiat-band host density as the hybridization 
function. In addition, the fiat-band case of the SIAM is very well understood P] and is 
thus helpful to explain features found in our model. We also take the host bandwidth 
to be much larger than any other bare energy scale and use a hybridization function 
which is constant, A(a;) = Aq. Our goal is then to compute the spectral density in the 
relevant energy window —W/2 < uo < W/2 with W/2 > U/2,Aq. For all numerical 
results presented here, the energy scale is set by Aq = I/tt. 

3. Method 

In this work, we employ the DMRG technique [S^ which is a reliable numerical method 
for one- dimensional systems. We use the standard DMRG method to calculate ground- 
state properties and the DDMRG method [^3] to calculate dynamical properties. In 
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Figure 1. One-dimensional lattice configuration for applying the DMRG algorithm 
to the two-impurity problem. The solid circles denote the impurity sites and the open 
circles denote the host band. 



order to carry out our calculations, we consider -|- 2 electrons in a system consisting 
of noninteracting bath sites {N even) and two impurity sites. The electron density 
is (nt) = (n^) = N/2 + 1. 

The (D)DMRG calculations can be performed on finite lattices only, i.e., we must 
discretize the host band and carry out (D)DMRG calculations for finite number N of 
host band eigenstates corresponding to energies Ek {k = 1, ...,N), and then extrapolate 
the results to a continuous host band (A^ oo) if needed. Choosing a discretization of 
the host band, i.e., selecting the N band state energies e^, should be done appropriately 
depending on what is to be obtained. 

The Hamiltonian is, however, somewhat unsuited for a DMRG treatment 
because it includes hopping terms that are long-range. For example, the system size 
is limited to < 60 for typical calculations even when several thousand density- 
matrix eigenstates are kept in the DMRG procedure, the maximum possible on current 
workstations. Therefore, we first transform the Hamiltonian into a linear chain with 
nearest-neighbor hopping only, 

i(T j(T ja 




The new fermion operators Cj^ correspond to electronic states in the host band and are 
related to the original representation by a canonical transformation 

Cj. = T.M,kL. (3) 

k 

The orthogonal matrix Mjk, the diagonal terms aj, and the nearest-neighbor hopping 
terms are calculated using the Lanczos algorithm for tridiagonalizing a symmetric 
matrix starting from the initial vector Mi^i^ = Vk/V with V"^ = J2k^k- -^^^ ^ 
hybridization function symmetric about u = 0, the diagonal terms aj vanish. The 
Hamiltonian Q describes two impurities coupled to one end of a one- dimensional chain 
representing the host band states (see Figure Q). This transformation enables us to 
handle a system with up to A^ ~ (9(200) bath states using the (D)DMRG method. 

Furthermore, for efficient treatment, we introduce even- {p = e) and odd- {p = a) 
parity impurity orbitals dpa = {di^ ± d2a-)/V^ as in Ref.l9. With Ri — R2 = 0, the 
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Hamiltonian (Q) is transformed to 

H =Ho + Hu (4) 
^0 = Vj2{dlco. + h.c.) + ^ A,(c]^c^.+i, + h.c.) 

+ t,,Y.idlA. - dlA.) (5) 

HU = ^[(^eT + ^oT - l)(^ei + rial - 1) 

+ {dl^do^ + dl^d^^){dlAi + dlidei)] , (6) 

where ne{o)a = dli^o)ade(o)(7 ■ Note that only the even-parity orbital hybridizes with the 
noninteracting bath states directly. When U = 0, the even- and odd-parity orbitals are 
completely separate. In our model, they will be mixed only via the Coulomb interaction. 

4. Results 

4.I. Static properties 

We begin our discussion with some static properties of the system, namely, the electron 
density at the impurities and the spin-spin correlation function between the impurities, 
calculated with the standard DMRG method. Here we apply a logarithmic discretization 
scheme of the host band, Ek = {W/2)A~'' (with A > 1 and k = 1, 2, ■ ■ ■ , A^/2), as 
usually used in Wilson's renormalization group method [2H|, because we are interested 
in the case of large host bandwidth as well as in having dense bath states around the 
chemical potential for quantitative accuracy. Typically, we use = 38 bath states with 
A = 1.5 and W = IOOttAq and keep m = 2000 density-matrix eigenstates in the DMRG 
procedure. In some cases, systems with up to A^ = 58 and W = IOOOttAq are used to 
extrapolate the results to ^ oo. 

4-. 1.1. Electron density at the impurity The average electron densities (he) and (fio) 
for even- and odd-parity orbitals of coupled impurities, respectively, are displayed in 
Fig. |2l^a) and (b) for different values of U as function of ti2. Since the two-impurity 
sites are equivalent, we have (fii) = (^2) = {{fie) + (^o))/2. Note that (ni) ((^2)) can 
take values between and 2 due to the charge degrees of freedom at the impurity sites, 
in contrast to the two- impurity Kondo model for which {fii) = {7x2) = 1 always. 

When ti2 = 0, {Uf.) = {ho) = 1 for all interaction strengths due to particle-hole 
symmetry. The local densities are drastically affected by tu. We find a discontinuous 
transition at a critical value t^ = ti2,c which is dependent on U. The critical value ti2,c 
is zero at f/ = and becomes larger with increasing U. In the large U limit, it saturates 
at ti2,c ~ ttAq. We show the critical value ti2,c as a function of t/ for — > cxd in the 
inset of Fig. Efc). 

At f/ = 0, {fio) drops to zero for infinitesimally small ti2(< 0). In this limit, there is 
no hybridization with the conduction band and no coupling to the even-parity level, i.e.. 
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Figure 2. Local density on the even- (a) and odd- (b) parity orbitals of tlie coupled 
impurities as a function of ti2 for several U values, (c) Spin-spin correlation function as 
a function of ti2 and U. Inset: the critical values of the transition, ti2,c, as a function 
of U. The dashed lines correspond to ti2,c = U/4: and ti2,c = 1-13 (see text). 



the odd-parity orbital forms a completely local state at the Fermi level ep iff ti2 = 0. An 
infinitesimally small ti2 shifts the state above ep, i.e., ti2,c = 0^. For small but finite U, 
the situation is similar and the transition is caused by the competition between tu and 
U. At ti2 = 0, the odd-parity orbital is split by the Coulomb interaction [see Eq.®] into 
two states with energy difference ~ t//2, the lower and the upper Hubbard band (we call 
them the LHB and the UHB, respectively) located at a; ~ ±f//4. A finite ti2 has the 
effect of shifting the LHB to higher energies by |ti2|. The odd-parity orbital is occupied 
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by one electron if the LHB is below ep and almost vacant if the LHB is above ep- This 
leads to ti2,c ~ U/i and a large discontinuity (~ 1) of {fio) at ti2,c- This anticipated 
behavior is in fact observed in our DMRG results. As U increases, the odd-parity 
orbital begins to hybridize indirectly with the host band via the interaction term (jH)). 
The discontinuity, therefore, becomes smaller and goes to zero a.s U —>■ oo. Moreover, 
we find that ti2,c has almost no [/-dependence in the large-f/ regime. This implies that 
the physics behind the transition for large U is different from the competition between 
ti2 and U at small U. Before we discuss a possible mechanism in the next paragraph, let 
us briefly discuss the behavior of the even-parity orbital. For all U values, {"he) increases 
slowly as a function of ti2 because the even-parity orbital hybridizes directly with the 
host band, thus experiencing a strong broadening. Hence, the discontinuity of (rie) at 
ti2 = ti2,c is smaller than that of (ho). We can also see the largest discontinuity, which 
means that the charge fluctuation on the even-parity orbital is largest, for intermediate 
U values (f/ ~ SttAq). We also note that (ho) = 1 independent of U below ti2,c, and 
(ho) — >■ and (he) ^2 as ti2 — >■ oo for finite U. 

4- 1-2. Spin-spin correlation between impurities We now investigate the spin-spin 
correlation between the two impurities, ^5*1 ■ 62^. The results are shown in Fig. |2fc) 
as a function of ti2 for different values of U. Before we discuss the results, let us first 
identify the different types of magnetic interactions present in our model. First, we have 
the c-f exchange interaction Jcf, which is the antiferromagnetic interaction between an 
electron on the impurity sites and conduction electrons. For a single impurity, Jcf leads 
to the Kondo effect and a local spin singlet as the ground state. Within the standard 
Schrieffer- Wolff mapping, the value for Jcf is given by Jcf = SAq /Ues, where Ues is 
the effective Coulomb interaction and the impurities are both occupied by one electron. 
The c-f exchange is effective only when an electron is localized on the impurity so that 
Ucs = U/2 and Jcf = 1QAq/U. Second, the conduction electrons mediate the RKKY 
interaction Jrkky- For our particular setup, we obtain a ferromagnetic interaction 
between two electrons on the impurity sites. As usual, Jrkky is obtained as second- 
order process, i.e., Jrkky ~ Jcf ~ C(l/f/^). Third, the model exhibits a direct exchange 
interaction Jex due to the coupling ti2, which is an antiferromagnetic interaction and is 
given by the standard expression J^x = 4t^2/^- 

We expect that for tu < ti2,c ferromagnetic correlations due to the RKKY 
interaction are dominant, i.e., (^Si ■ 6*2^ > 0. On the other hand, when ti2 > 
ti2,c, antiferromagnetic correlations due to the exchange interaction are stronger, i.e., 
(^Si • 5*2^ < 0. Such a transition from ferromagnetic to antiferromagnetic correlations at 
ti2,c is in fact found at all interaction strengths, as can be seen in Fig.|2fc). The absolute 
value of (^Si ■ 82^ increases with increasing U and reaches the maximum possible value 
as ?7 ^ 00, which means that one electron is localized on each impurity in the U 00 
limit. 

Let us now consider the [/-dependence of the critical value ti2,c- The transition 
is driven by the charge degrees of freedom in the weak-coupling regime. However, the 
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spin degrees of freedom play an essential role in the strong-coupling regime. Taking 
into account that the spin-spin correlations between the impurity sites change from 
ferromagnetic to antiferromagnetic at ti2,c, we may expect that the competition between 
the RKKY and direct exchange interactions is the origin of the transition. If we 
take Jex = </rkky as the criterion for the occurrence of the transition, we obtain 
ti2,c ~ 0{1/VU). Thus, ti2,c would go to zero as [/ cxd, which is obviously inconsistent 
with the DMRG results which show almost constant ti2,c as a function of U. 

Up to now, we have not taken into account Jcf, which leads to a competition 
between the formation of local Kondo and non-local singlets, offering a quite different 
mechanism for the transition. In this case, the criterion to obtain ti2,c is Jcf = Jc^, 
which gives ti2,c = \/4/7r ~ 1.13. This result is indeed consistent with our findings. 
Thus, the transition in the strong coupling regime can be interpreted as competition 
between local singlet formation due to the Kondo effect and non-local singlet formation 
due to the direct exchange introduced by ti2. Note, however, that, in contrast to the 
general folklore, the boundary is not set by Tx^Jcf) = Jex, but by the direct comparison 
of the bare energy scales. We also find that the total spin is 5" = 1 for ti2 < ti2,c and 
S* = for ti2 > ti2,c • 



4-2. Dynamical properties 

In this section, we study the spectral density for impurities in the TIAM. The impurity 
one-particle Green function for even- and odd-parity orbitals can be written as 

\ H - Eq + uj -irj I \ Eq - H + UJ + 17] / 
{rj — > 0+), where Eq is the ground-state energy and (...) represents a ground-state 
expectation value. The impurity spectral density for each parity is then obtained as 

Dp^{(jj) = -- sgn{uj)lmGp„{uj) = Ap^{uj) + Bp„{uj) (8) 

71 

with 

Ap„{uj < 0) = hm (rft dr,a) 

Bp^{uj > 0) = hm ( i . dl\ 

and Aprj{uj > 0) = Bp^^u < 0) = 0. The spectral density fulfills the sum rule 

Dp^{uj)du = 1. (9) 



Note that the spectral densities for both impurities are the same, i.e., Di„{uj) = D2a{^^), 
and are equal to [Doo-(i^) + -Deo-(t^)]/2. 

The standard DMRG algorithm [311 HO] can be used to calculate the ground-state 
properties as shown in the last section. In particular, the ground-state wave function 
l^'o) and the ground-state energy Eq can readily be obtained. To compute dynamic 
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properties such as the impurity Green's function ((7j) we use the DDMRG • This 
approach is based on a variational principle. One can easily show that for rj > and 
fixed frequency u, the minimum of the functional 



(^Eo + u - hY + f]^' 



d. 



pa 



^0 



with respect to all quantum states |\[') is 



-rj 



pa 



Eo + uj-H) +ri 



-dl 



pa 



(10) 



'111 



The functional minimum is related to the convolution of the spectral density with a 
Lorentz distribution of width r] by 

Wi^^;^) = --nvB'Uuj). (12) 



if one substitutes d^^^ for dj,^, 



A similar result is obtained for the spectral density (j9j) 
—uj for oj and A^^ioj) for B^„{ui) in the above equations. 

The DDMRG method consists essentially of minimizing the functional (jlUj) 
numerically using the standard DMRG algorithm. Thus the DDMRG provides the 



spectral densities A^^ 



pa 



.00 



and B;^{u) 



for a finite broadening rj. The full spectral 



density (jS)) convolved with the Lorentz distribution 



D 



pa 



V 



7r[(cj 



UJ 



/\2 



(13) 



is given by the sum of A^^ 



pa 



UJ 



and B^i^iuj) 



The real part of the Green's function 
can be calculated with no additional computational cost but is generally less accurate. 
The necessary broadening of spectral functions in DDMRG calculations is actually very 
useful for studying continuous spectra or for doing a finite-size scaling analysis fHK] . 

What we really would like to obtain is the spectral density in the rj = limit. This 
can be done by carrying out a deconvolution of the DDMRG data [36]. In theory, a 
deconvolution amounts to solving (fT^ for Dpa-{uj) using the DDMRG data on the left- 
hand side. We also know that the broadened spectral density of the impurity system 
on an infinite lattice (A^ —>■ oo) is usually almost identical to the spectral density of 
the discretized impurity system (A^ < oo) if rj > Ae. Therefore, one can make the 
approximation that the DDMRG data for D^^{u) describes the broadened spectral 
density for N oo and can then solve ()13|) approximately under the condition that 

n 

D. 



pa 



.UJ 



pa 



UJ 



is the exact spectral density of the TIAM. For instance, one can require that 
is a continuous and relatively smooth function. To obtain quantitatively accurate 
spectra after deconvolution, we need to take r/ smaller than the width of the spectra in 
ri = 0. 
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Figure 3. Spectral density at [/ = and W = 207rAo. The open circles denote De„{u}) 
calculated with a constant host band discretization for N = 59, Ae « 0.347rAo, and 
77 = O.SttAo, then deconvolved. Solid circles denote Deai^) calculated with variable 
discretization Ae w O.OGTvrAo and r/ = O-IttAq around the peak. Solid and dashed lines 
are the exact solutions (|14|l with broadening 77 — O.IttAq and (|15|l without broadening, 
respectively. 



4-2.1. noninteracting case The spectral density at f/ = can be calculated exactly and 
provides a good test to demonstrate the accuracy of our method. The exact spectral 
density for the odd-parity orbital is a a (5-function at u; = —tu 



while for the even-parity orbital we obtain a Lorentzian of width 2Aq centered at uj = ti2 



In Fig. El we show the spectral density calculated with the DDMRG for U = 0. On the 
scale of Fig. El there is no visible difference between our numerical results and the exact 
results. The deconvolution technique is not useful for obtaining a divergent function 
such as a ^-function, so we introduce a finite broadening 77 = CIttAq into Eq.lfT^ 
and compare it to our 'bare' DDMRG spectra convolved with same rj. Note that the 
local Coulomb interaction U is always treated numerically exactly in the density-matrix 
renormalization and thus does not affect the accuracy of the method directly. 

4-2.2. weak- coupling regime In the weak-coupling regime {U <^ 47rAo), the physical 
properties are still similar to those of the noninteracting case. The spectral density 
of the impurities calculated with the DDMRG for ti2/Ao = 0, O.IStt, and 0.257r at 
U = ttAq is shown in Fig. |31 The critical coupling ti2,c is 0.2057rAo. Let us first look 
at the spectrum for the even-parity orbital, Deo-('-t'). At tu = [see Fig. jH^a)], Dea{uj) 




(14) 




(15) 
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Figure 4. Spectral density of tlie coupled impurities for the odd-parity orbital Doa{^) 
and for the even-parity orbital Df.a{u)) with (a) ti2 = 0, (b) ti2 = O.ISttAo, and (c) 
ti2 = 0.257rAo at /7 = ttAq and W = 207rAo. Dashed lines denote Decr{'-^) calculated 
with a constant host band discretization for N = 58, Ae w 0.347rAo, and ry — O.SttAo, 
then deconvolved. Solid lines denote Doa{'-^) calculated with a variable discretization 
for N = 118 [N = 70), O-OItt < Ae/Ao < 0.457r (0.00337r < Ae/Aq < 1.987r) and a 
constant broadening r/ = 0.027rAo (77 = COOSttAq) for ti2 = and CISttAq (0.257rAo), 
then deconvolved. Insets: expanded view around the Fermi level lj — Q. 



is basically a Lorentzian of width ~ 2Ao centered at = 0, but there appear small 
shoulders around u ~ ±0.87rAo(~ U) and uj ~ ±1.87rAo(~ 2U) due to the Coulomb 
interaction When ti2 increases, the central peak is shifted towards lower energies 
by |ti2| while maintaining its shape [see Fig. El^b) and (c)]. This is consistent with the 
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gradual increase of the local density on the even-parity orbital (ng) as a function of t^ 
and also with the small discontinuity at ti2,c, as shown in Fig. Efa). We next turn to 
the spectrum for the odd-parity orbital, Docr{i^)- As long as f/ = 0, this orbital has no 
direct hybridization with the host band, so that Doo-(co') consists of a localized state, 
i.e., a single 5-function peak. When a small f/(< ttAq) is introduced at t^ = 0, this 
peak splits into two peaks located around ~ ±f//4 due to the effective repulsion on the 
odd-parity orbital Ues = U/2 [see Eq.®]. They correspond to the LHB and the UHB. 
In other words, the odd-parity orbital is half-filled and a Mott-Hubbard gap opens. The 
peaks are still very sharp but are no longer exact ^-functions, as we can see in Fig. ^a), 
because the odd-parity orbital couples 'indirectly' to the host band via the even-parity 
orbital. As ti2 increases, the two peaks are shifted towards higher energies by |ti2|, but 
their separation remains ~ U/2. The LHB and UHB become sharper (broader) while 
retaining their respective weights. In addition, no spectral weight is transferred to the 
gap. When the LHB reaches the Fermi level ep {u = 0) at ti2 = ti2,c(~ U/A), the 
transition occurs. When t^ > t^^c [see Fig. IHc)], only one peak is present, and it is 
above uj = and is very sharp. Actually, there must be some spectral weight below 
uj = because (ho) is not exactly zero [see Fig. Efb)]. We therefore find that the odd- 
parity orbital behaves as a nearly localized state for all ti2 and that the charge degrees 
of freedom in the odd-parity orbital play a crucial role for the transition. 

4-2.3. strong- coupling regime We now consider the spectral density in the strong- 
coupling regime {U ^ 47rAo), which corresponds to the so-called Kondo regime in 
the SIAM. For the TIAM, however, the situation is more complex due to competing 
interactions, i.e., the RKKY interaction, the Kondo (or c-f exchange) effect, and the 
exchange interaction between impurity sites. In Fig. we show the spectral density 
calculated with the DDMRG for ti2 = 0, O.GttAo, and ttAq at U = ISttAq. The critical 
inter-impurity hopping ti2,c ~ O.SttAq is slightly smaller than that obtained in Sec. II 
because the calculations were done for a finite host bandwidth here. Let us first look at 
the even-parity spectral density Deaiyj). Below ti2,c, we can see a sharp peak at = 
in Deaijjj), which satisfies the Friedel sum rule D^aiyJ = 0) = l/(27rAo). This means 
that the conduction electrons form a spin-singlet (Kondo) state with electrons in the 
even-parity orbital. The width of this peak at a; = becomes smaller exponentially with 
increasing U . We can thus state that the properties of Df-aiyj) below ti2,c are similar to 
that of the SIAM which is characterized by the Abrikosov-Suhl resonance at = and 
the Hubbard satellites around uo ±[//2. Moreover, we notice that the physics of the 
TIAM is quite different from that of the two-impurity Kondo model, where no Kondo 
effect is observed j^, at least for the case of equivalent impurities. Note also that in 
the TIAM, states with only one electron on the two impurity sites can still have a large 
weight (~ 3% for U = ISttAq) in the eigenvector of the ground-state even in the Kondo 
regime, in contrast to the two-impurity Kondo model. 

When we increase ti2, the shape of De„{uj) is hardly changed and only a weak 
transfer of spectral weight from above ep to below ep occurs, consistent with the 
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Figure 5. Spectral density of tlie coupled impurities for the odd-parity orbital Doa{^) 
and for the even-parity orbital Dea-{'-^) with (a) ti2 = 0, (b) ti2 ~ O.GvrAo, and (c) 
ti2 = ttAo at U = ISttAq and W — 407rAo. Dashed lines denote Z?eo-(<^) calculated 
with a variable discretization for N = 94, Ae sa O.OOISttAo, and r] = 0.0027rAo, then 
deconvolved. Solid lines denote Doa{'^) calculated with a variable discretization for 
N = 94, 0.00687r < Ae/Ap < 1.987r and a constant broadening -q = O.OIttAo, then 
deconvolved. Insets: expanded view around the Fermi level uj = Q. 



behavior of {fie)- As long as ti2 is smaller than ti2,c, the quasi-particle peak stays pinned 
at cij = and maintains its height (l/(7rAo)), while the Hubbard satellites stay located at 
\uj\ <U /2 with width > 2Ao. Above ti2,c, on the other hand, the Kondo peak vanishes. 
The local spin is now screened due to the formation of a non-local singlet between the 
impurities due to the dominant exchange interaction, i.e., the scattering channels leading 
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to the Abrikosov-Suhl resonance are not active any more. Simultaneously, the 'effective' 
hybridization between the host band and the even-parity orbital becomes weaker so that 
the width of Hubbard satellites becomes narrower. Our numerical results indicate that 
the Hubbard satellites turn into Lorenzians with width 2Aq and weight 1/2 located at 
\uj\ = U/2 in the limit of ti2 — > oo and U — >■ oo. 

Next, let us turn to the odd-parity spectral density Doa-{^^)- At ti2 = 0, we find four 
peaks. The two prominent peaks located at ^ ±7.5(= ±7rAof//2) can be identified 
with the Hubbard bands of the electrons localized in this orbital. Since no direct 
hybridization to the band states exists, these peaks are sharper than the corresponding 
ones in Dea{oj), i.e., the broadening is introduced indirectly via the Coulomb interaction. 
Rather more interesting is the appearance of a structure at the Fermi energy consisting of 
two narrow peaks separated by a gap A ~ 0.36 ~ Jc/ = 16/ (vrt/) ^ 0.34 for ti2 = 0. As 
ti2 increases, this feature and, in particular, the gap prevails, although its size decreases 
(see insets to Fig. EI). Note, however, that the gap edges remain symmetric with respect 
to the Fermi energy, while the spectral weight is larger below the Fermi level as long as 
ti2 < ti2,c- For ti2 > ti2,c5 a rearrangement of spectral weight from below to above the 
Fermi energy takes place, while the size of the gap does not change noticeably. 

We interpret this structure as a replica of the Kondo resonance induced indirectly 
by the interactions between the even and odd channel. However, adding or removing 
electrons in the odd-parity orbital would at least break one "Kondo bond" in the even- 
parity channel, i.e., cost an energy ~ Jcf, explaining the appearance and size of the gap 
for ti2 = 0. Note that, in contrast to the weak-coupling regime, this explanation connects 
the gap to the spin rather than to the charge degrees of freedom. As ti2 increases, the 
gap becomes smaller because the direct exchange interaction Jex between the impurity 
sites competes with the Kondo effect, effectively reducing the local moment, which is 
screened by the band states, and hence the corresponding spin gap. However, we find 
A > Jcf — Jcx, which would be the naive expectation, because the formation of dynamic 
spin correlations between the impurities due to Jex introduces a further contribution to 
A. Thus, the gap will be always finite and develops a minimum close to t^.c, where 
•Jcf = Jex- For ti2 > ti2,c5 the impurities form a non-local spin singlet and the spin 
gap will scale with Jex, i.e., increase again slowly with increasing tu- An estimate for 
ti2 > ti2,c yields A ~ Jex = 4t^2/^ ~ 0.26, which is in rough agreement with the 
numerical value A ^ 0.16 obtained from Fig. \^c). 

4.2.4- intermediate- coupling regime Finally, we present the spectral function for 
U ~ 47rAo, i.e., in the intermediate-coupling regime. Note that from the point of view 
of the SIAM this value already resides within the "strong-coupling" regime delimited 
by = 2 (the effective hybridization for the even-parity channel is 2 Aq) . In Fig. IHl 
we show the spectral density calculated with the DDMRG for tu = 0, O.SttAq, and 
ttAq at [/ = SttAq. The transition here occurs at ti2,c ~ 0.657rAo. As already 
mentioned above, the even-parity spectral density D^aiuj) below ti2,c resembles that of 
the SIAM. The central peak at a; = starts to become narrower and its spectral weight 
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Figure 6. Spectral density of tlie coupled impurities for the odd-parity orbital Doa{^) 
and for the even-parity orbital Deai^) with (a) ti2 = 0, (b) ti2 = O.SttAo, and (c) 
ti2 = ttAq at U = SttAq and W — 207rAo. Dashed hnes denote Deai^^) calculated 
with a constant discretization for N = 58, Ae « 0.347rAo, and t] ~ O.SttAo, then 
deconvolved. Solid lines denote Doa{^) calculated with a variable discretization for 
TV = 118, O.OItt < Ae/Ao < 1.457r and a constant broadening r/ = 0.027rAo, then 
deconvolved. Insets: expanded view around the Fermi level lo — 0. 



is increasingly transferred to the high-energy range with increasing U. However, the 
three-peak structure typical of the strong-coupling regime is not yet fully developed; 
the Hubbard bands appear only as visible but shallow shoulders around u ~ ±U/2. 
Moreover, in contrast to the SIAM, the central peak here does not fully reach the 
Friedel limit, its height being slightly lower than I/tt. This behavior was also seen in 
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previous NRG studies [22] and in the mean-field approach j2Sl , and is connected to the 
non-local magnetic correlations induced by the RKKY exchange, which is stronger for 
smaller U. As in the strong-coupling regime, Dea{'^) hardly changes with increasing ti2 
until ti2,c is reached, where a dramatic redistribution of spectral weight connected to the 
formation of a non-local singlet due to the direct exchange introduced by ti2 appears. In 
particular, as in the strong-coupling regime, the Kondo peak has vanished completely. 

The behavior of the odd-parity channel in Fig. IHl is also similar to the strong- 
coupling limit. We again find a gap in the spectrum, which remains symmetric about 
= for all ti2, and observe a similar but much more pronounced change in the 
distribution of spectral weight, as in Fig. The gap, too, initially decreases until ti2,c 
is reached and then increases again. We believe that the physics behind this behavior is 
essentially the same as in the U oo limit, although the energy scales associated with 
the spin gaps induced by the different exchange mechanism in particular now cannot be 
written down explicitly. However, we expect them to be larger than in the limit of large 
U, which is indeed what we observe in Fig. IHl 

5. Conclusion 

In this paper, we have presented static and dynamical properties of the two-impurity 
Anderson model at half filling for vanishing impurity separation. In this limit, the 
otherwise rather complex model becomes considerably simpler, nevertheless retaining 
most of its interesting physical aspects. In particular, the competition between different 
types of magnetic correlations such as Kondo, RKKY and superexchange is preserved. 
We employ the (dynamical) density-matrix renormalization method to calculate static 
properties (filling, magnetic correlations) and one-particle spectra. For the latter, we 
obtain results that are convoluted with a Lorenzian, which we, however, can deconvolute 
with good accuracy, as demonstrated for the exactly solvable case U = 0. 

While this model can, in fact, easily be studied by e.g., Wilson's NRG, an evident 
advantage of the DDMRG is that one obtains an accurate description of spectral features 
on all energy scales at the cost of the resolution of exponentially small structures. In the 
NRG, on the other hand, exponentially small scales can be readily resolved, however, 
at the expense of accuracy at intermediate and high energies. Since we are interested 
here in describing both the features emerging at intermediate (~ Jrkky, Jex) and high 
energy scales (Hubbard bands ~ U/2) as well as at a possible small Kondo scale, we 
feel that the DDMRG is more useful for the present study. 

In the weak-coupling regime, U < 27rAo, we observe a transition between a 
situation with weakly ferromagnetically coupled impurity spins to a situation with weak 
antiferromagnetic correlations as function of inter- impurity hopping tu, as is apparent 
from the behavior of (5^i ■ S2) in Fig. Efc). From the spectral functions discussed in 
Fig. I3J we can furthermore infer that this transition is primarily driven by the charge 
degrees of freedom in the odd-parity channel of the two-impurity system. For U > 27rAo, 
on the other hand, the transition is into a state with rather strong antiferromagnetic 
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correlations. At the same time, the behavior of the spectral functions, in particular in the 
odd-parity channel, changes considerably. Although the transition is still accompanied 
by a change in occupancy in the odd-parity channel, this change is visibly reduced. 
Furthermore, the gap in the odd-parity spectrum is always pinned to the Fermi energy 
and is of the size of the typical magnetic exchange interactions, which points to the 
spin degrees of freedom as the driving force of the transition. Finally, the existence 
of a Kondo peak in the even-parity spectrum for ti2 < t^^c, which abruptly vanishes 
for ti2 > ti2,c, must be taken as evidence that the transition is in fact driven by the 
competition between Kondo and direct exchange in this regime. It is quite important 
to note that the actual transition occurs at the point where Jc/ = Jex and not when 
TxiJcf) = -'cx, and also that the energy scales appearing in the spectra in Figs. Eland 
ini are in fact related to these "intermediate" quantities rather than the actually much 
smaller Kondo scale. 

Of course, the current investigation is restricted to a special limit of the TIAM, 
namely, that of vanishing impurity separation. The results for this case clearly show that 
a more thorough investigation of this model is still necessary. In particular, we believe 
that such an investigation must carefully study the relation between the different energy 
scales inherent to the problem. The method to study the TIAM should be chosen so that 
it can at least resolve accurately intermediate, i.e., of the order of the effective exchange 
interactions, and low energy scales, i.e., of the order of the Kondo scale. Evidently, we 
cannot expect that any method can handle both regimes equally well when the energy 
scale of the latter is exponentially small, but we believe that the DDMRG can at least 
treat the intermediate situation when Tr- <^ J^x but is still resolvable with the method. 
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